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We analyze solutions to Friedmann-Robertson-Walker cosmologies in Brans- 
OS ■ Dicke theory, where a scalar field is coupled to gravity. Matter is modelled by a 

7-law perfect fluid, including false-vacuum energy as a special case. Through 
a change of variables, we reduce the field equations from fourth order to 
second order, and they become equivalent to a two-dimensional dynamical 
system. We then analyze the entire solution space of this dynamical system, 
and find that many qualitative features of these cosmologies can be gleaned, 
including standard non-inflationary or extended inflationary expansion, but 
also including bifurcations of stable or unstable expansion or contraction, 
noninflationary vacuum-energy dominated models, and several varieties of 
"coasting," "bouncing," "hesitating," and "vacillating" universes. It is shown 
that inflationary dogma, which states that a universe with curvature and 
dominated by inflationary matter will always approach a corresponding flat- 
space solution at late times, does not hold in general for the scalar-tensor 
theory, but rather that the occurence of inflation depends upon the initial 
energy of the scalar field relative to the expansion rate. In the case of flat 
space (k = 0), the dynamical system formalism generates some previously 
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known exact power-law solutions. 
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I. INTRODUCTION 



It was first recognized by Guth that inflation, i.e., a period of superluminal expansion 
in the early universe, solves many of the problems of the standard model of cosmology, 
a.k.a. the "hot big bang" theory These include the smoothness/horizon problem, the 
flatness/age problem, the structure problem and the massive relics problem. (A review of 
these problems, and their resolution in an inflationary scenario, is given in [Q.) However, 
in the original scenario proposed by Guth, one finds that the deSitter expansion induced 
by vacuum energy density (or, equivalently, by a cosmological constant) in ordinary general 
relativity is impossible to exit from via barrier penetration in a way that permits sufficient 
inflation to solve the aforementioned problems, while still reheating the universe afterwards 
||. In particular, it can be shown that if the tunneling transition rate for the inflaton 
field is small enough to allow for sufficient inflation, then it is too small to nucleate true- 
vacuum bubbles fast enough to collide and reheat the universe, and the phase transition to 
a radiation-dominated era is never completed. A number of possible solutions have been 
proposed for this "graceful exit problem," including an extremely flat inflation potential 
("new inflation"), so that inflation occurs as the field "rolls down" towards its true- vacuum 
state, without any barrier penetration being necessary ||]; and including "chaotic inflation," 
in which only certain regions of the universe inflate to large size before the true vacuum is 
restored ||. 

One such proposal, "extended inflation" ||, based on the Brans-Dicke (BD) theory of 
gravity |7|], solves the graceful exit problem while still allowing the phase transition that 
ends inflation to be completed via bubble nucleation, by employing power-law rather than 
exponential expansion. This permits eventual complete reheating even if the tunneling 
transition rate for the inflaton is small enough to ensure that inflation continues at least 
until the relevant cosmological problems are solved. 

In its original form, this proposal fails due to the creation of unacceptably large microwave 
background perturbations by collisions between big bubbles |J. To avoid this "big bubble 
problem" , the coupling constant must take on a value in the range u < 25 during inflation, 
in conflict with limits set by recent experiments of u > 500 |J. This conflict may be 
circumvented by further alterations of the gravitational sector, such as a potential for the 
BD scalar [jlTJ. In addition to the big bubble problem, however, one finds that the spectrum 
of density perturbations produced by quantum fluctuations of the BD scalar during inflation 
is distorted away from scale-invariance, so that COBE measurements of the spectral index of 
the perturbations put a further constraint on the allowed range of us during inflation. Taken 
together, these constraints appear to rule out many extended inflation models ||11|| . Models 



which survive include those with a variable coupling constant such as hyperextended inflation 



12| , and hybrid models which include both a first-order phase transition and a period of 



slow-roll |[3] . 

Recent interest in BD theory is also motivated by its relation to low-energy bosonic string 
theory. We can write the action for BD theory as 



Sbd = J (Tx^g ^-4>R + ^^-^ + 167r£ m J , (1.1) 
where R is the scalar curvature, u is the dimensionless BD coupling constant, and <ft is the 
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BD scalar field. The low-energy effective action for bosonic string theory can be written in 



the form 14 



S ef f = J d A x^e~ 2 " (r + 4$'"$„ - ^H 2 ^j . (1.2) 

Here R is again the scalar curvature, but now $ is the dilaton field and H 2 = H llvp H lll/p , 
where H^ vp is the totally antisymmetric three-form field. If we make the substitution = 
e -2 * we see that, considering only the scalar and tensor components and neglecting the 
matter, the two theories are identical in the low-energy limit if u — — 1. It should be 
emphasized, however, that the two theories differ in their couplings of the scalar field to 
the other matter. The renewed interest in scalar-tensor gravitation has led to several recent 
investigations into the generation of exact solutions for cosmology in such theories [30], as 
well as to some qualitative studies of the models which result f3"T|,|32|j . 



There is considerable interest, therefore, both from inflationary cosmology and from 
string theory, to revisit cosmology in the scalar-tensor gravity of Brans and Dicke. Many 
papers were of course written on BD cosmology in the 1960's and 1970's ||T5|-19|; interest 



then declined, partly because observational tests of gravity put increasingly stringent lower 
limits on the BD u in the present universe, and partly because fundamental scalar fields fell 
out of theoretical favor. However, string theory restored the theoretical credibility of the 
scalar field, because of the necessary appearance of the dilaton, leading to the realization 
that the scalar field can, and indeed must, gain a mass and decouple from present-day long- 
range physics. So, one's point of view on BD cosmology has changed: it does not apply at 
the present time, but may apply during eras of inflation, or even of quantum gravity, in the 
early universe. Differences from general relativity may be relevant to the viability of string 
theory at energies below the Planck scale, and to inflationary scenarios which incorporate 
scalar fields coupled to the other matter. 

However, if one is to use the scalar-tensor theory in the place of general relativity in 
early cosmology, a natural and necessary first step is the delineation of the full spectrum of 
cosmological models which are possible in the theory. With the addition of an adjustable 
coupling constant and a scalar field requiring two initial conditions, it should come as no 
surprise that a much richer variety of qualitatively distinct models are available here than in 
general relativity. The purpose of this paper is to consolidate these models into a complete 
and coherent representation of their character. To our knowledge there has been no previous 
attempt to do this. An essential practical tool here is a new transformation of variables that 
reduces the system of equations from fourth order to second order, greatly simplifying the 
analysis. 

After this work was completed, we learned from D. Wands that there is substantial 



overlap between our work and his Ph.D. thesis P3| , including the reduction of the equations 
to a second order system. 

The paper is organized as follows: after writing down the equations for homogeneous 
and isotropic cosmology in the scalar-tensor theory in Sec. II, we reduce those equations to 
a 2-dimensional dynamical system in Sec. III. The analysis of the dynamical system, which 
contains descriptions of all of the qualitatively distinct models allowed by the theory, is 
performed in Sec. IV. In Sec. V we use the dynamical system formalism to rederive some 
previously known exact power-law solutions to the theory in flat space, and Sec. VI presents 
conclusions. 
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II. COSMOLOGY IN BD THEORY 



The field equations for BD theory can be found by varying the action, Eq. (|1 . 1|) , with 
respect to the metric and the scalar field. They are 



uj- 



+ R = 



(2.1) 



and 



57T< 



-lj. , 
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;9»u<i>' p <l>,p) 



+ 



g^n<f>), (2.2) 



where 







-gdx. 



Taking the trace of Eq. ( |2.2|) and combining the result with Eq. (2.1), we find 



S7T 



3 + 2a; 



TV 



(2-3) 



(2.4) 



Eqs. Q2.2| , [2.4|) may be viewed as the fundamental equations for BD theory. We take as the 
background metric the homogeneous and isotropic Friedmann- Robertson- Walker (FRW) line 
element, which may be written in co moving coordinates as 



ds 2 = -dt 2 + a 2 (t) 



dr 2 



1 — kr' 



+ r 2 dn 2 2 



(2-5) 



where a(t) is the scale factor and k = 0, ±1 is chosen to represent the universe as a space of 
constant zero, positive or negative curvature, respectively. Then we find that the nontrivial 
components of the field equations, Eqs. (|2.2j , P-4|) , are 



+ 



+ 



1 d <l 
cr dt 



k _ /2cu + 3 
^ ~ 12 
8tt 



3 + 2o; 



"3^ 



(2.6) 
(2.7) 



where we have written only the time component of the metric field equation, since the 
spatial components are dependent on it through the Bianchi identities. Assuming a perfect 
fluid form for the stress-energy tensor, i.e., T^ u = diag(p,p,p,p), we also find the usual 
conservation equation (the time component of T^ v '. u = 0): 



P = -3-(p + /»). 



(2.8) 



Eqs. ( |2.6| - |2.8D are the equations for FRW cosmology in BD theory. Note that we have a 
system of one second-order and two first-order differential equations, so that we need four 
pieces of initial data to completely specify a solution in this theory, as well as an equation 
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of state p(p) and values of u and k. By contrast, when solving for the Friedmann models in 
general relativity, we have two first-order equations and hence need only two pieces of initial 
data. The difference stems from the fact that in BD theory, <p~ l acts like a time-varying 
gravitational constant obeying the second-order Eq. (|2.7|) , so that we need two additional 
initial conditions in order to close the system. An interesting consequence of this is that 
even after specifying initial conditions for a, (f) and p, we are left with a one-parameter 
family of solutions parametrized by the initial value of the rate of change of the scalar 



field. Indeed, as has been noted by other authors |15| , p0| , the character of the solutions at 
early times depends solely upon the value of the constant a 3 <p at the initial singularity. At 
later times, however, most solutions become indistinguishable from those where a 3 <fi = at 
the singularity. Accordingly, there has been a concentration in the literature on the three- 
parameter family of solutions characterized by a 3 (O)0(O) = 0. However, one must bear in 
mind that there has been a loss of generality when this boundary condition is imposed. 

We turn now to theoretical considerations regarding the possible values of the coupling 
constant uo. First, inspection of Eq. ( |2.6| ) reveals that at early times, the curvature term 
vanishes relative to the other terms provided only that a < as a — > 0, i.e., provided that 
the initial epoch in the history of the universe was not inflationary. If we make the further, 
seemingly reasonable assumptions that p > and <fi > 0, then we see that we must have 
u > —3/2 in order to satisfy Eq. ( |2.6| ) for all times. Perhaps a more compelling argument 
for this constraint on uj arises when BD theory is reformulated in a way that preserves the 
form of the Einstein equations (the so-called "Einstein frame"), by treating the scalar as a 
dynamical field which couples to all of the other matter When this is done, one finds 
that the stress-energy for the scalar field is 

rjf = J^(<^ - (2-9) 

so that assuming homogeneity, 

From this, we see that the energy density of the scalar will be nonnegative in the Einstein 
frame only for uj > —3/2; thus, for example, one can only prove the existence of black holes 
or of a cosmological singularity for u in this range. As a result of these considerations, we 
take to > —3/2 in what follows. 



III. REDUCTION OF THE FIELD EQUATIONS TO A DYNAMICAL SYSTEM 

In this section, we change variables in order to transform the fourth-order system specified 
by Eqs. ( |2.6| - 12.8|) into a pair of coupled second-order equations in which, however, only first 



derivatives of the new variables appear. The resulting equations may then be analyzed 
using the established tools of dynamical systems theory. First we switch to conformal time 
dr = dt/a(t), and define the new variables 

(3.1) 




6 



a 



2u; + 3\ 6' = </>' 
12 ~~ = " _ ' 



where primes represent derivatives with respect to conformal time, and 



.4 
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Next, we parametrize the equation of state by writing 

P = (7 - !)P- 



(3.2) 



(3.3) 



(3.4) 



Since the three-form field is necessarily present in the low energy limit of string theory, it is 
of interest to consider how it might bear on the present study. If we assume a homogeneous 
and isotropic universe, then it follows that 



Hq^ — 0, 

#123 = Kt). 



If we then use the fact that 



T, 



a/3 



-2 6(y/=gC) 

a/3 



-9 Sg 



(3.5) 



(3.6) 



to calculate the stress-energy for the three-form under these assumptions, we find that it 
can be modelled phenomenologically as a perfect fluid with equation of state p — p. Then 
we have, for example, 
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0, false- vacuum energy; 

1, pressureless dust; 
4/3, radiation; 

2, 3-form field. 



(3.7) 



Now we can rewrite Eqs. (|2.6| - |2.8| ) as 



(3 2 + k 



(a 2 6')> 



P = -37-P 

a 



iirpa 2 a 2 
36 + ~A 



(1 - 3 7 /4) 



3A 



a 



(3.8) 

(3.9) 
(3.10) 



Note that although false-vacuum energy may play a role in the dynamics of the theory, 
the cosmological constant A does not explicitly appear in the field equations as we have 
written them. Unlike in general relativity, where the inclusion of A is exactly equivalent to 
the presence of vacuum energy, in BD theory the two are qualitatively distinct. This is the 
case because whereas in GR the vacuum energy enters the field equations through the term 
8irGp v (which one can equate to A), in BD theory it enters through the term 87r p v j '6, clearly 
not in general a constant. Others have considered the case A / with no matter present 
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22|; in a separate paper, the present authors will consider the system with both A ^ and 



additional matter present [23[ . 

Combining Eq. ( |3.8| ) with Eq. ( |3.9| ) and using the definitions of (3 and a, we can derive 
the expression 

a' = (1 - 3 7 /4)03 2 + k - a 2 /A) - 2/3(7. (3.11) 



Next we take the conformal time derivative of Eq. ( |3.8| ), and use this in conjunction with 
Eq. ( EUOl) and Eq. ( pi]) to obtain 

P' = (1 - 3 7 /2)(/3 2 + fc) - (3 - 3 7 /2)(x7A (3.12) 

Eqs. ( |3.11| , |3.12j ) constitute a dynamical system (in conformal time) for the variables (3 and 
a, and as such it is conducive to both quantitative and qualitative analysis. 



IV. QUALITATIVE ANALYSIS 

Although it is certainly agreeable to have some exact solutions to the system specified by 
Eqs. ( p.6| - p.8|) , such solutions typically correspond to a particular choice of the initial data. 



Furthermore, it is difficult to synthesize the conglomeration of known particular solutions to 
the theory into a picture which reveals its overall character. It is therefore desirable to obtain 
a more complete understanding of the solution space, which may be accomplished through a 
qualitative analysis of the dynamical system ( |3.11| , |3~T2j ) in the j3-a plane, accompanied by a 
"translation" of the results back into the physically meaningful a-<p plane. Before beginning 
the analysis, however, we shall make some general comments regarding the system. 

First, let us consider what an equilibrium point of the system represents. Whenever 
(3' = a' = 0, this implies that a = (3 — a/2A is a constant, so that 

«(*) = (ft - ^) t + ao. (4.1) 
Also it is straightforward to show that at such a point, 

0(t) = CA [((3 - a/2A)t + aoY /[A{ ^ a/2A)] (4.2) 
unless j3 = &/2A, in which case 

<f>(t) = o exp( ( xt/Aao). (4.3) 

Hence we see that a fixed point in the (3- a system represents either a solution where the 
scale factor changes linearly with time while the scalar varies as a power of time, or a static 
universe with exponentially varying scalar. 

Next, it will be of interest, as we consider the various fixed points of the system, to 
determine whether or not a given point lies in the physical domain p > 0. We can check 
this using Eq. ( |3.8| ), which reads 



( 3* + k = a*/A+ 6:Tpa2 



S 



so that (assuming > 0,) any point (/?o,co) satisfying 



Po 2 + k > a 2 /A (4.4) 

will lie in p > 0, whereas points satisfying 

(3 2 + k< a 2 /A (4.5) 

lie in p < and thus do not represent physical solutions. 

A brief discussion regarding equations of state is in order at this point. Of particular 
interest is the question of what equations of state will lead to inflation in this theory. As was 
mentioned in Section II, BD theory can be reformulated through a conformal transformation 
so as to preserve the form of Einstein's equations, while viewing the scalar field simply as a 
part of the matter (albeit a peculiar part) which couples dynamically to every other field, 
and hence to the particle masses. Using this form of the theory, we can write the usual 
equations of general relativity for homogeneous, isotropic cosmology, i.e., 



ay k 8nGp 



,2 



- = — =— (Ptot + 3p tot , (4.7) 
a 6 

where the subscripts are a reminder that these quantities include contributions from the 
scalar field in addition to the "regular" matter. As we saw previously [cf. Eq. (|2.10 )1, such 
contributions are given by 

We see that inflation will occur only when p tot < —ptot/3- For models which are not dom- 
inated by the scalar field, this just reduces to the usual condition for inflation in general 
relativistic cosmology, p m < —pm/3 (7 < 2/3). It is also clear, however, that the con- 
tribution of the scalar field to the stress-energy tends to act against inflation (and, more 
generally, against expansion), so that inflation may not occur if the scalar field dynamics 
become sufficiently dominant, regardless of the equation of state of the other matter. 

Finally, the applicability of the singularity theorems of Hawking and Penrose to 
BD theory is of interest. In general relativity, these theorems prove the existence of a cos- 
mological singularity for the Friedmann- Robert son- Walker models, provided that Einstein's 
equations and the strong energy condition are all satisfied. For a perfect fluid, this latter 
condition reduces to p + 3p > 0, i.e., only for inflationary matter is there the possibility 
of a cosmological model without a singularity. To apply the theorems to BD theory, we 
must transform to the Einstein frame, where the Einstein equations are formally preserved. 



Under this transformation, one finds that El 



a(t) - Gv 1/2 r 1/2 a(t), (4-9) 

from which it follows that 
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BD 



a -) • (4.10) 

Oj) e 



Hence we see that an expanding universe in the Einstein frame is anything with (in our 
notation) (3 > in the BD frame. As a result of this, a model with ordinary matter 
which appears to contract and then smoothly reexpand in the BD frame will nevertheless 
be consistent with the singularity theorems so long as (3 > during the transition, for then 
the model is in fact always expanding in the Einstein frame and the singularity in the past 
would be due to the infinite value of the scalar. 

Now let us proceed with the analysis. Rewriting the results of section III, we start with 

P'=(l- 3 i)(F + k)-(3-^)r-), (3.12) 




° =\ l - -r)\P + k - — ) - (gjjD 

where, as before, 

(3.1) 

V 12 ' ■ — - (m 

and we are in conformal time dr = dt/a. It is convenient to divide the equilibrium points 
(Po, Co) into two categories: 

(i) (3 2 + k = 0, a = (k = 0, -1) 

(it) /?o = ±iRl-3 7 /4), a„ = ± t p(l-3 7 /2) (Vfc) (4.11) 




where 



H = 3A(1 - 7 /2)(l - 3 7 /2) - (1 - 3 7 /4) 2 . (4.12) 



First consider = 0, in which case all of the equilibrium points, regardless of the values 
of 7 and lu, take on the values (/3 ,o"o) = (0,0). According to our recent discussion, this 
point represents the static solution a = do, <p = 0o- Inspection of the field equations 
also reveals that this solution is only valid in the true vacuum p = 0. However, although 
the solution appears at first glance to be unphysical, it may also be viewed as the usual, 
asymptotically approached solution d — >■ 0, a — > oo for a flat, expanding universe. The 
fact that the rate of change of the scalar field approaches zero at large times is simply the 
statement that the strength of the gravitational interaction becomes nearly constant, as it 
must if the theory is to agree with what we observe in the universe today P3| . Fig. la shows 



the numerically integrated solution curves in the (3-a plane for a flat, matter-dominated 
(p=0) universe, where we have arbitrarily selected A — 1 (u — 9/2) for the plot. The 
shaded areas represent unphysical regions of negative energy density, and the dark line 
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represents the transition between expansion and contraction, this transition being forbidden 
for positive energy densities in the present case. The solutions to the right of (3 = are 
initially expanding, and we see that they asymptotically approach the equilibrium point 
(/3o, 0"o ) = (0,0), an interpretation of which has already been given. The solutions to the 
left of /3 — are initially contracting, and will continue to do so until a singularity is 
reached. We see that the solutions approach and depart from equilibrium along a particular 
path, which represents a known power-law solution to the flat-space field equations. We 
shall consider such exact solutions, and their relation to the general class of cosmological 
models, in Section V. For the time being, we simply mention that the attractive nature of 
such solutions near equilibrium is a general feature of the theory, and that these solutions 
typically correspond to the particular choice of initial data a 3 = at the singularity. 

A qualitative change occurs when the coupling constant takes on a value in the range 
—3/2 < uj < 0. Fig. lb is similar to Fig. la, except that here we take uj = —1/2. Note 
that now the line representing a = exists in the physical regime where p > 0. In fact, one 
sees that models exist where the universe is initially contracting, but then makes a smooth 
transition to expansion without meeting a singularity. This behavior has been noted by 



other authors |19| , and models with this feature are sometimes called "bounce" universes. 
It should be mentioned that such models are present for oj — — 1, the value corresponding 
to the low-energy limit of string theory. 

Next we examine the flat space models with an inflationary equation of state. Fig. lc 
shows the solution curves for k = and p = —p, where again we arbitrarily select uj = 
9/2. The attracting nature of the particular solution represented by the line (3 = 3a is 
quite apparent, and this is in fact the inflationary power-law solution which has been found 



previously by other authors ||26|| . Here all expanding universes asymptotically approach this 
state, and once again all contracting universes must reach a singularity, although we see that 
as they do so, they approach the true-vacuum solutions (represented by the dark lines at 
the border of the shaded region). These are examples of what may be termed '^-dominated' 
models; i.e., models where the scalar field dynamics overwhelm the effect of the other 
matter and thus attract the universe towards true-vacuum behavior. Nevertheless, these 
models must contract until a singularity is reached. 

Again the character of the solutions changes in the range —3/2 < u < 0, in a manner 
similar to the noninflationary case. Here, however, not only do we see that there exist 
"bounce" models, but we also find that inflation is no longer possible in this range of uj. 
Rather, all models approach the equilibrium solution a = 0, a — > oo as they would for 
noninflationary equations of state; another example of '0-domination'. This behavior is 
illustrated in Fig. Id. In the range < uj < 1/2, inflation is still not possible, although the 
"bounce" models are absent. 

Now let us consider the case of the open universe, k — — 1. First consider the equilib- 
rium points in category (i), i.e., (/3 ,cr ) = (±1,0). Here the corresponding solutions are 
a(t) = ao±t, <f)(t) = 0o which are, again by inspection of the field equations, valid only for 
p = 0. The physical interpretation of the point (/3 , o"o) = (+1, 0) is straightforward; we 
can view it either as representing a particular vacuum solution, or as the asymptotically 
approached endpoint of an open, expanding universe. The contracting solution represented 
by the point (/?o, o"o) = ( — 1, 0) does not lend itself to interpretation quite so easily; to aid us 
in our efforts we can perform a stability analysis. This proceeds by linearization, i.e., we 
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check the eigenvalues of the Jacobian matrix for the system ft' = f(ft,a)] a' = g(ft,a). In 
cases where the eigenvalues of the Jacobian all have nonvanishing real part, the fixed point 
is called hyperbolic and we can determine its stability from the signs of those real parts: if 
the real part of each of the eigenvalues is negative at a given equilibrium point, the solution 
is stable at that point; if the real part of each eigenvalue is positive, or if the real part of 
one eigenvalue is positive and that of the other is negative, then the solution is unstable at 
that point. Finally, if the real part of any of the eigenvalues is zero at a point, then the 
point is called nonhyperbolic and its stability in the neighborhood of that point cannot be 
determined by this method P7|. Hence we write 



2/3 (l-3 7 /2) -T( 3 " 
2ft {l - 37/4) - 2a -^(1-37/4) 




+ ... , (4.13) 



Jacobian 



where £ = ft — fto and r\ = a — a , and diagonalize the Jacobian case by case. In this 
manner we find that for the equilibrium points in the k = — 1 models under consideration, 

if (ft Q , cr ) = I I , then 



A : =2(l-3 7 /2), A 2 = -2 
A! =2(3 7 /2-l), A 2 = 2 



( stable, 7 > 2/3 
[ unstable, 7 7^ 2/3 



(4.14) 



As expected, the point (1, 0) represents a stable solution for noninflationary, open uni- 
verses. The point (—1, 0), although an equilibrium point of the dynamical system, in general 
represents an unstable and thus unphysical contracting solution. Hence a small perturba- 
tion will send the universe along a more rapidly contracting trajectory as shown in Fig. 2a. 
In any case, all models with uj > 0, negative curvature and noninflationary equations of 
state must, if initially expanding, approach equilibrium at d = 0, a — > 00 or, if initially 
contracting, reach a singularity. If the equation of state is inflationary, i.e., if 7 < 2/3, then 
the expanding equilibrium solution represented by (1,0) becomes unstable, and in the range 
w > 1/2 we do in fact see inflation for these models. Although the existence and stability 
of the equilibrium points (±1,0) for k = — 1 do not depend upon the value of uj, we again 
find models which can contract and smoothly reexpand if —3/2 < uj < 0. This behavior 
is illustrated in Fig. 2b for models dominated by pressureless dust. For k — — 1 models 
with —3/2 < uj < 1/2 and inflationary equations of state, inflation does not occur despite 
the fact that the equlibrium solution at (1, 0) is unstable. This is illustrated in Fig. 2c and 
Fig. 2d. In order to better understand these models, we must move on to an analysis of the 
remaining class of equilibrium points. 

Let us consider, then, case (ii) of Eq. (|4.11|) , i.e., the equilibrium points with o~o 7^ 0. 
Referring back to that equation, we find that such points will exist only if the condition 
k/^Ujj) > is met, which implies that 



/2u + 3\ 1 / 2 



< 7 < - 
' 3 



1 + 



f2u + 3\ 1 / 2 



iuj 



(4.15) 
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if k = — 1, and 



7 > 



1 + 



[2uo + 3\ 1/2 



or 7 < 



1 - 



2u + 3\ 1/2 



(4.16) 



if k — +1. We can thus think of the curves 



7 



1 ± 



/2lu + 3\^ 2 



as bifurcation curves in the parameter space (uj, 7). By this one means simply that there are 
equilibrium points of the dynamical system which move in or out of existence as those curves 
are crossed. However, using Eq. ( p.8|) and Eq. (|4.11|) , we find that these fixed points will lie 
in the domain p > only if 7 < 2/3, i.e., only for matter with an inflationary equation of 
state. Fig. 3 shows the bifurcation curves in the u—j plane. The region between the curves 
contains new equilibrium points for k — — 1; the regions exterior to the curves contains such 
points for k — +1. The unshaded region represents that part of the parameter space with 
equations of state which can support the solutions represented by these new fixed points and 
still maintain p > 0. It is evident that the only bifurcation of physical interest occurs for 
7 = 0, i.e., with vacuum energy dominant. In that case, one finds that two new equilibrium 
points appear at uj = 1/2. If the curvature is positive (k = +1), then these equilibrium 
points exist in the regime uj > 1/2; if the curvature is negative (k = —1), then they exist 
in the regime uj < 1/2. Let us first complete our analysis of the negative curvature models. 
Taking k = -1 in Eq. ( jOTD 



we have 



(A, °v) = (±(1 - 3A)- 1 / 2 , ±A(1 - 3A)-^) (uj < 1/2), 



for which we find from the Jacobian that 



Ai < 0, A 2 < 
X 1 > 0, A 2 > 



(4.17) 



(4.18) 



Here we see that one of the equilibrium points is stable (a sink) while the other is unstable 
(a source), as we have already seen explicitly in Fig. 2c and Fig. 2d. Let us consider the 
nature of the stable solution. Using the definitions of (3 and a, one finds that the solution 
at this point is given by 



a(t) = a + (£/2)t; 
0(t) = Ca 2 , 



(4.19) 



where 



E = (1 - 3A) 



-1/2 



(4.20) 



and, by assumption, p = p vac = const. Hence, although not inflationary, the universe will 
expand forever as we would expect for a model with negative curvature. Also, one finds that 
in order to satisfy the field equations, the value of p vac must be given by 
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ZAC (2uj + 3)C 

Pmc ~ 4tt(3A - 1) 47r(l-2a;)' 1 J 

where C is the integration constant in the the solution for <ft, which we can write as C = 
0o/oq. Hence we see that for this solution, the vacuum energy is completely determined by 
the initial conditions of the scale factor and the scalar field, and by uj. 

Finally let us consider the models with positive curvature, i.e., those with k = +1. The 
analysis of the equilibrium points is, in light of our previous discussion, particularly simple. 
There are no equilibrium points in category (i) of Eq. (|4.11|) , and the equilibrium points 
from category (ii) which lie in the domain p > exist only for inflationary equations of state 
and only for uj > 1/2. Then using Eq. (|4.11|) , we find 

(A), *o) = (±(3A - I)" 1 / 2 , ±A(3A - l)- 1 / 2 ) ( W > 1/2), (4.22) 

and diagonalization of the Jacobian at those points yields 

i.e., both of the equilibrium points are unstable. So for models with positive curvature, 
we have no physically interesting equilibrium points regardless of the equation of state and 
the value of uj. Nevertheless, we may determine the character of the solutions from their 
evolution in the /3—a plane. In general, one expects that such models should recollapse 
unless inflation occurs. Indeed, one finds that for noninflationary equations of state, i.e., 
7 > 2/3, and in the range uj > 0, this is in fact the case. However, as before, models with 
— 3/2 < uj < and those with 7 < 2/3 have features unique to this theory. To begin with, 
one finds that in the range —3/2 < uj < 0, there exist noninflationary models which do 
not collapse, despite having positive curvature. This is illustrated for pressureless dust in 
Fig. 4a. Note that the initial conditions must be carefully tuned for such models, and in 
particular that o < always, so that is a monotonically decreasing function of time. 

For inflationary matter and positive curvature, consider first the range —3/2 < uj < 0. 
Here one finds that inflation never occurs, although there exist models with carefully chosen 
initial conditions which expand forever. All other models, either initially expanding or 
initially contracting, collapse to a singularity. These models are illustrated in Fig. 4b. 
Next, we consider vacuum energy-dominated models with k — +1 but with < uj < 
1/2. Here we find that inflation never occurs, and furthermore that all models, regardless 
of initial conditions, collapse to a singularity. Finally, consider inflationary models with 
positive curvature and in the range uj > 1/2, where we have seen that there exist 2 unstable 
equilibrium points in the (3-a plane. Here we find solutions with several interesting features. 
First, in the range 1/2 < uj < 9/2, there are models which start out expanding, and then 
recollapse rather than inflate, including some which "coast" at a = before beginning their 
contraction. If uj > 9/2, then in addition we find that there exist "bounce" models which 
start out contracting, but then smoothly reexpand and inflate — these are similiar to general 
relativistic de Sitter spacetime in k = +1 coordinates. Also in this range of uj, there are 
models which are initially contracting, then smoothly reexpand, and finally recontract until 
a singularity is reached — we call these "vacillating universes." Indeed, there exist models 
which are initially expanding, then recontract, and at last smoothly reexpand and inflate. 
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H uj > 9/2, then in general one finds that the universe will inflate if and only if its initial 
conditions lie in the region bounded by the stable manifolds of the two equilibrium points. 
On the other hand, any model which has initial conditions which lie between either of the 
stable manifolds and the true vacuum in the j3-a plane will collapse back to a singularity. The 
fact that the stable manifolds extend arbitrarily far back in the negative /3-direction gives 
rise to yet another feature unique to this theory: an arbitrarily rapidly contracting universe 
will slow down, smoothly reexpand and inflate if the energy of the scalar field is chosen just 
right. These models represent a particularly glaring form of "0-domination" : whether the 
universe collapses to a singularity or inflates depends solely upon the energy of the scalar 
relative to the expansion rate. Fig. 4c illustrates this with uj = 10.5. The division of the (3-a 
plane into inflationary and noninflationary cells by the stable manifolds is apparent. Fig. 4d 
is a closeup of the neighborhood of the equilibrium point in the expanding regime of Fig. 4c, 
and one sees explicitly the "bounce" and the "vacillating" universe models. The limiting 
cases of "vascillation" are clearly present as well. For example, a universe which slows down 
in its expansion (or contraction), and then speeds back up may be termed a "hesitation" 
model, whereas one which just stops but then continues its expansion after a finite time may 
be called a "coasting" model. The potentially drastic effects of the scalar field can also be 
seen in the special case of u — 9/2, where one finds that a necessary condition for inflation 
to occur is P > cr, which implies that 

(a) > 2(0) 

(However, even if this relation is satisfied, the universe may recollapse if 0/0 is negative 
enough initially.) For this value of a;, a universe which is initially static will begin to expand, 
following one of the stable manifolds towards its unstable equilibrium point, although a small 
perturbation in the metric or the scalar will cause either collapse to a singularity or inflation, 
depending upon the direction of the perturbation. 



(4.24) 



V. EXACT SOLUTIONS FOR FLAT SPACE 



Many exact solutions to Eqs. ( pT6| - p.8[ ) have been found for the case of flat space (k = 0), 



some of which yield power-laws for the evolution of the scale factor, the scalar field, and the 
energy density PS|| . Here we find exact power-law solutions which are parametrized by the 
equation of state, and show that these include all possible power-law solutions to the field 
equations. 
If we define 

A-§, (5.1) 

then using Eqs. ( fj. 1 1| , PTT2|) with k = 0, we can derive the expression 

dX (A 2 -A)[3A(l- 7 /2)-(l-3 7 /4)] 



d\n(3 A(l - 3 7 /2) -3A 2 (1- 7/2) 



(5.2) 



Note that any constant value of A which is a root of the cubic polynomial in the numerator 
of the right-hand-side will satisfy this equation. Those roots are 
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A = ±VJ (5.3) 
(1-37/4) 

" 3(W2)' ( } 

the first of which represent the true vacuum solutions, and the last of which represents a 
solution with matter (or vacuum energy) present. Before finding these solutions explicitly, 
let us remark that the initial condition (pa 3 — > as a — > is implicit in most of the solutions 
that we can generate with this formalism. To see this, first note that if follows from the 
definitions of (3 and a that if A is a constant, then 



(5.5) 



where we have set <fr = <fto when a = ao, and 

T- A/2 



(5-6) 



Then we can write 

0a 3 ~ da Q+2 , (5.7) 
from which we see that if a ~ t Cl at early times, then 0a 3 — > as a — > so long as 

d > (5.8) 
3 + a 

One can easily verify that this relation is satisfied for all physically reasonable perfect fluid 
equations of state, i.e., for 7 = 0, 1, or 4/3, unless: 

- 5/6 < uj < -1/2 for 7 = 0; 

-4/3 < uj < -1 for 7 = 1. (5.9) 

To see that the A = constant solutions give us all possible power-law solutions, assume 
a power-law solution <fi oc t p , a oc t q . It then follows from the definitions of (3 and a that 

i.e., A = constant. 



A. True- Vacuum Solutions 

If we set A = ±v^4, we obtain the true- vacuum solutions. This can be seen by using 
Eq. ( |5.5|) and setting p = and k = in Eqs. (|2.6| , |2.7| ), upon which we find that for 
A = ±y/A, Eq. (|2.6|) is satisfied identically, whereas Eq. ( |2.7| ) can be integrated to yield the 
power-law solutions 
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a(t) = a [l + (tAo)] (±2 ^- 1)/(±6v/3 - 1) , 

0(t) = o [l + (tAo)] 2/(±6 ^ 1) (5.11) 

where we have set (p(t = 0) = 0o and a(t = 0) = ao, and to is arbitrary. These solutions are 
not valid for uo = —4/3, since the exponents corresponding to the positive roots in Eq. ( |5.11| ) 
blow up there. The special true- vacuum solution corresponding to uj = —4/3 was found by 
O'Hanlon and Tupper in [|I^] to be deSitter in nature, and may be written in the form 

a(t) = a exp(t/t ), 

0(t) = o exp(-3t/t o ). (5.12) 



B. Solutions with Matter Present 



Now consider the solutions where p ^ 0. Then from T 0v = we have 

P o h 

- = -37- 
p a 

( a \ - 3 t 
P = Po[ — 



(5.13) 



where we have set p = po when a = a . Using Eqs ( |5.5| , |5.13j ) in Eqs (|2.6| , |2.7|) , we find that 
they become 



[l + a/2) 2 - Ao? 



a\ 2 8irpo(a/ao) 



-37— a 



30o 



and 



+ 



a 87r(4 — 37)^0(^/00) 



-37 



(5.14) 



(5.15) 



a 3 + 2uj 

Again setting a{t = 0) = a and <p(t = 0) = (f> , Eq. (|5.14j) immediately integrates to 

a(t) =a [l + X t} 2/{a+3l) (5.16) 

where 

\ 1/2 



X 



a + 37 



57rp 



2 J V30 o [(l + "/2) 2 - Aq 2 [ 
Now using Eq. (|5.5|) and Eq. (|5.16|) , we see that 

Also, using Eq. (|5.13|) , we find 

p(t) = Po (i+xr 67/(a+37) 



(5.17) 



(5.18) 



(5.19) 
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Eqs. ( |5.16| - |5.19|) specify all power-law solutions to flat space FRW/BD cosmology. 



Now we can use these results for the various types of matter. First consider vacuum 
energy, where p = —p vac is a constant, 7 = 0, A = 1/3, and a = A/{2uj + 1). Then using 
Eqs. (|5TI^5U9|) , we find 



a(t) = oo(l + X tT +1/2 

( j ) (t) = M^ + Xt) 2 (5.20) 



with 



y 2 = ^lt Pvac 

Xvac o (6cu + 5)(2cu + 3) 1 ' 

in agreement with the known power-law solutions |j,26|. Note that if uo < 1/2, then a < 
and inflation does not occur. This is noteworthy, as we see that the bare low-energy limit of 
string theory, corresponding as it does to a value of uj = —1, does not permit an inflationary 
solution for flat space. In Section IV, we saw that this result also holds for a universe with 
curvature. 

Next consider a MD era, where p = 0, 7 = 1, A = 1/6, and a — + 1). We find 

a(t)=a (l + X tf" +2)/{3 " +4) 

p(t)=p (l + X t)- 3 ^ +2)/ ^ +4) (5.22) 

where 

2 47rp f3uj + 4 



= ^7 toJ (5 - 23) 

again in agreement with the known power- law solutions f?|,pO|j. 

Finally, for a RD era, we have p = p/3, 7 = 4/3, A = 0, and a = 0. From this we find 

a(t) = a (l + xtf 2 
0(t) = 0o(l + Xt)" 1/2 , 

pit) = p (l + XtV 2 (5.24) 



where 



&» = (5-25) 
300 

Note that in this case, the scale factor and energy density evolve just as they do in ordinary 
general relativity, independent of the coupling constant to. This is the case because the 
scalar field has been effectively decoupled from the other matter due to the vanishing trace 
of the stress-energy tensor [cf. Eq. (|2.4j )1. 
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VI. CONCLUSIONS 



The resurgence of interest in scalar-tensor gravity is justified, in part, by the formal sim- 
ilarity of the theory to low-energy string theory, and by the power-law inflationary solution 
which ameliorates the graceful exit problem of some general relativistic models of inflation. 
However, in order to create a viable cosmological model with the scalar-tensor theory, one 
must understand, to as full an extent as possible, what models are available in the theory and 
how those with qualitative differences may be distinguished from one another. Here we have 
shown that there exist a plethora of models which differ substantively from those present 
in the standard general relativistic cosmology. In general, the models are parametrized by 
initial conditions of the scalar field, the value of the BD coupling constant, an initial value 
of the Hubble parameter, and an initial density and equation of state of the matter present. 
The first two parameters are not, of course, present in general relativity. 

For flat space, in addition to the standard models which start with a singularity at the 
big bang and then asymptotically approach a = (or the time-reversal of this behavior) as in 
general relativity, we have found "bounce" models, which pass smoothly from contraction to 
expansion. Such models are present for all types of matter in the range of coupling constant 
uj < 0. In addition, we have found that a flat-space inflationary solution does not exist in 
this theory in the range uj < 1/2, irrespective of initial conditions. 

For a universe with negative curvature, we again find a set of models which have an 
analogue in general relativity, in this case those which start with a singularity and then 
asymptotically approach the solution a ~ t (again the time-reversal of this exists with 
t — ► —t). And again, we find "bounce" models when u < 0, as well as the absence of an 
inflationary solution for uj < 1/2. We find that in this latter range of uj, all expanding neg- 
ative curvature models with inflationary-type matter present approach a particular solution 
at late times, where a(t) ~ t and <ft(t) ~ t 2 . 

For a universe with positive curvature, the standard models are once more present, but 
anomalous behavior occurs both in the range —3/2 < uj < and for inflationary matter. 
Again inflation will not occur if uj < 1/2, in which case all models with vacuum energy 
dominant collapse to a singularity only in the range < uj < 1/2. If uj is negative, there 
exist models which perpetually expand for all types of matter, including vacuum energy. 
If vacuum energy is present and uj > 1/2, we find several classes of models which have 



no analogues in general relativity without the addition of a cosmological constant p9 
First, inflation still does not occur for certain ranges of initial conditions, those ranges 
being determined by the locations of the stable manifolds of the dynamical system in our 
formalism, which in turn are determined by the value of uj. Next we find "hesitating" and 
"coasting" models, which start with a big bang and either slow down and "hesitate" before 
continuing their expansion, or "coast" at a quasistatic value of a(t). When uj > 9/2, we find 
models which expand from a big bang, contract for a time and then reexpand ( "vacillating 
universes" — the time-reversed solution also exists) and models which pass from "deflation" 
to inflation, i.e., universes which, although contracting arbitrarily fast, turn around and 
ultimately inflate. 

Of particular interest are those models with nonzero curvature that are dominated by 
vacuum energy which, contrary to the well-accepted inflationary dogma, do not approach a 
corresponding flat-space solution at late times. Indeed, we have seen that there exist, for all 
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values of the coupling constant, vacuum energy-dominated models with curvature for which 
there are no corresponding flat-space solutions at all. In such cases, the initial energy of the 
scalar field relative to the expansion rate sometimes plays a key role in determining whether 
inflation will occur, a fact which appears to have been largely overlooked in the literature. 
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FIGURES 



FIG. I. The evolution of solutions in the (3-a plane for flat space, (a) Matter is 7 = 1 
(pressureless dust); and uj = 9/2, representative of uj > 0. Solutions to the right of the line 
P — a/2 represent expanding universes; those to the left are contracting universes. The 
shaded regions require negative energy density and so are disallowed physically, (b) Same 
as (a), except uj = —1/2. The line ^ = has moved into the regime of positive energy 
density, so that some models pass smoothly from contraction to expansion, a general feature 
of models with —3/2 < uj < 0. (c) Flat space 7 = (vacuum energy dominated) models, 
with iv = 9/2. All initially expanding universes become inflationary — they flow out along 
the prominent channel to the left — regardless of initial conditions when k = and iv > 1/2. 
(d) Flat space, 7 = (vacuum energy dominated) models with iv = —1/2, representative of 
— 3/2 < uj < 0. Inflation does not occur, although all initially expanding universes continue 
to expand forever. 

FIG. 2. Some models with k = — 1 (negative curvature), (a) 7 = 1 (pressureless dust) 
with uj > 0. (b) 7 = 1 with iv < 0. (c) 7 = (vacuum energy) with uj < 0. (d) a closeup of 
the stable equilibrium point in (c). In (a) uj = 9/2; in (b-d) iv = —1/2. 

FIG. 3. Bifurcation diagram in the parameter space tv-j. The shaded region is disallowed 
physically as it requires p < 0. The regions exterior to the curves contain bifurcated equi- 
librium points for k — +1; the region interior to the curves contains bifurcated equilibrium 
points for k = —1. The bifurcation for k — +1 and 7 = (vacuum energy) is physically 
permissible and occurs when iv = 1/2. 

FIG. 4. Some models with k — +1 (positive curvature), (a) 7 = 1 (pressureless dust) 
and uj = —1/2, representative of —3/2 < uj < 0. Some initially expanding models do not 
recollapse, despite positive curvature, (b) 7 = (vacuum energy dominant); and uj = —1/2, 
representative of —3/2 < uj < 0. Most expanding solutions recollapse; a few expand forever, 
but inflation never occurs, (c) 7 = (vacuum energy dominant); and uj = 10.5, represen- 
tative of uj > 9/2. Most expanding universes inflate, but some do not. Here an arbitrarily 
rapidly contracting universe will slow down, stop and inflate if the energy of the scalar rela- 
tive to the expansion rate is carefully chosen, (d) A closeup of an unstable equilibrium point 
in (c). The stable manifold separates models which inflate from those which recollapse. Also 
apparent are "vacillation" universes which start out contracting, then briefly expand, and 
finally recontract to a singularity. 
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